Insulator, conductor and commensurability: a topological approach 
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I discuss a topological relation of the conduction property of a many-particle system on a periodic 
lattice at zero temperature to the energy spectrum. When the particle number per unit cell is an 
irreducible fraction p/q, an insulator must have q low- lying states of energy o(l/L) in one dimension 
and of energy o(l) in two dimensions, where L is the linear system size. 



An important question in the condensed matter 
physics is that under what conditions the system be- 
comes a conductor. While we do not know, of course, the 
complete answer, there is a characteristic of the system 
which may be naturally related to the conduction. In a 
many-particle system on the periodic lattice, the known 
insulating states such as band and Mott insulators can be 
understood as a "locking" of the particles by the lattice. 
It requires the particle density to be commensurate with 
the lattice, namely the filling factor (the number of par- 
ticles per unit cell of the lattice) to be an integer. When 
the filling factor is not integral but is a (simple) ratio- 
nal number, the particles may still be locked by spon- 
taneously breaking the translation symmetry and thus 
effectively enlarging the unit cell. When a small number 
of particles are "doped" into the system, they become 
mobile carriers, and we expect the system to be gapless 
and conducting. While this argument looks natural, it is 
far from trivial in the strongly interacting systems. It is 
needed to clarify how far this physical expectation holds, 
and what are the precise consequences. 

In the recent several years, it has been pointed out that 
the commensurability (in the sense of the particle number 
per unit cell) is related to the topological nature of the 
system. Extending the ideas by Lieb-Schultz-Mattis [lj 
(LSM) and by Laughlin Q, several statements on the 
energy spectrum related to the commensurability are de- 
rived j3, H ■ On the other hand, the naive expectation 
of the spontaneous breaking of the translation symme- 
try apparently contradicts with the "Resonating- Valence- 
Bond" (RVB) type gapped quantum liquid states. They 
have 1/2 particle per unit cell, but there is no long-range 
order which can be defined in terms of the local operator. 
(For recent examples, see e.g. Refs. H, El 0) However, 
this is still consistent H E| with the topological argu- 
ment as there is a "topological degeneracy" of the ground 
states under the periodic boundary condition. Moreover, 
the same underlying mechanism was used in a topologi- 
cal proof [TlLll2| of the Luttinger's theorem 0]. It was 
also applied to the Fermi sea volume of a Kondo lattice, 
which had been unamenable to rigorous analyses. 

On the other hand, there still remain several funda- 
mental questions, including (but not limited to) the fol- 
lowings. While the intuitive argument would relate the 



commensurability to the conduction property, the rela- 
tion has not been discussed in the topological context. 
Moreover, argument in Ref. for systems in more than 
one dimension relies on a rather nontrivial assumption 
that a finite gap does not collapse due to an Aharonov- 
Bohm (AB) flux 0. It is still necessary to clarify 
whether (or when) that assumption can be justified. In 
the present Letter, I shall address these questions by dis- 
cussing the Drude weight at zero temperature in view of 
the topological argument. 

The classification of the insulator and conductor has 
been discussed by several authors. Drude weight repre- 
sents the "free acceleration" of the current by the electric 
field, and thus gives a criterion to define the (perfect) 
conductor. The Drude weight D is defined by 
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where <j(u>) = <j(uj, q = 0) is the zero momentum compo- 
nent of the complex AC conductivity, a reg {uj) is assumed 
to be finite at lu ~ 0, and 5 is a positive infinitesimal. 
Kohn related the dependence of the ground state energy 
on the AB flux to the Drude weight [14| . His argument 
has been further developed by several authors, in par- 
ticular by Scalapino, White and Zhang |l5|. In addi- 
tion, Resta and Sorella £u5:| put forward an alternative 
but deeply related criterion of an insulating state which 
I will comment later. 

Since the analyses on the Drude weight in Refs. 0,0] 
and the topological argument on the commensurability 
are both associated with the insertion of the AB flux, one 
may expect they can be related to each other. In order to 
establish a concrete connection, let me first reformulate 
the Drude weight in the real time. I consider a system de- 
fined on a finite ci-dimensional periodic lattice of volume 
V = L X X Ly X . . . x Ld with periodic boundary con- 
ditions, where Lj denote the lengths in each directions. 
The lengths and the volume are measured in terms of 
the unit cell of the lattice. Let the filling factor (particle 
number per unit cell of the lattice) be v = p/q. I choose 
L x to be an integral multiple of q, and all other Lj's to 
be mutually prime with q. In taking the thermodynamic 
limit Lj —> oo eventually, v is kept constant. 

Suppose, for the time t < 0, no electric field is applied 
and the system was in the ground state. It can be chosen 
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so that it is an eigenstate of the x component of the 
total momentum operator P x . For < t < T a small 
constant electric field E is applied uniformly throughout 
the system. I then define the "time-dependent Drudc 
weight" 



D(t) = 



Et 



(2) 



where j(t) represents the current in a;-direction. Since 
the Drude weight represents the free acceleration of the 
current by the electric field, D(t) converges to D in the 
limit of t — > oo. Actually it can be proved, provided that 
the linear response theory holds. For a later convenience, 
I further define a weighted average D of D(t) over the 
time domain < t < T by 



(3) 



D = — I dt tD(t). 



In the limit of large T, D is dominated by the long-time 
asymptotic of D(t), and thus also should converge to D. 

One still has to be careful in using the present formu- 
lation, however. For a given lattice model and a given 
(finite) electric field E, (j(t)) cannot keep the t-linear 
behavior indefinitely. In particular, in a lattice model, 
{j(t)} would not exceed some finite upper limit. This is a 
manifestation of the non- linear effects. Generally, when 
(j(t)) becomes large, non-linear effects are important and 
the linear response theory breaks down. In order to apply 
the linear response theory (so that eqs. J2J and © gives 
the standard Drude weight), the current (j(t)) should be 
sufficiently small. 

Let us represent the electric field by the time- 
dependent vector potential A(t), which is uniform 
throughout the system, in the x direction. It is given 
by A(t) = Et. In fact, this A(t) is related to the AB 
flux $(t) which pierces through the "hole" of the torus. 
$(<) = L x A(t). The current j(t) is given by 
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where TC(&) is the Hamiltonian with the AB flux Since 
$ is the function of time, Hamiltonian is also a function 
of time. 



dH _ d<S> dH($) _ (M$) 
dt ~ dt (9$ x d<f> ' 
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Combining eqs. Q and l|Sj. we obtain the interesting 
relation j(t) = (l/EV)dH/dt. Now the averaged Drudc 
weight D is given by 



D 



T 2 E 2 V 
<5>(T) 2 V 



l dU \rli 



{{n) 



t=T 



(H) 



t=o 



(6) 



In the Schrodinger picture, (H)t=T should be understood 
as (V(T)\H($(T))\V(T)) where |*(t)) is the state vector 
at time t. This formula is exact and does not rely on the 
adiabatic theorem. 

If one takes T to infinity and $(T) to infinitesimal, the 
time evolution may become adiabatic (however, see dis- 
cussions below). In this case, the above formula reduces 
to the well known Kohn formula 



D = D 
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where Eq(Q) is the ground state energy as a function of 
the AB flux ($). 

Scalapino, White and Zhang pointed out a sub- 
tlety in the Kohn formula in particular for a gapless sys- 
tem in d > 1. Namely, many level crossing occur even 
for a very small AB flux $ in a large system, and the 
definition of d 2 E(Q)/d& 2 becomes actually rather sub- 
tle. They offer some definitions which would give the 
correct Drude weight. An alternative solution proposed 
here is to use the real-time formulation JJjJ instead of the 
Kohn formula Q . I fix the time interval T (and thus the 
"sweeping rate" of $) and take the thermodynamic limit 
Lj — > oo first. After taking the thermodynamic limit, I 
let T — > oo. There could be more and more level cross- 
ings and the time evolution is not necessarily adiabatic, 
as the system size is increased with a fixed T. With huge 
numbers of crossings in a large system, the final state 
|\I/(T)} may actually be very complicated 01 but D is 
still well-defined. By taking T — > oo limit after taking 
the thermodynamic limit, D should converge to the stan- 
dard Drude weight D. This procedure in the real time 
corresponds to, in the frequency space, taking uj — > 
limit after taking the thermodynamic limit with q = 0. 

Now let me choose $(T) to be the unit flux quantum 
$o = 27r/e 2 where e is the charge of the particle and I set 
% = c = 1. I note that, this choice corresponds to apply- 
ing a small electric field E = 0(1/ (TL X )). Even at t = T, 
the current is j(T) - DET = 0(l/L x ). Thus, it can be 
arbitrarily small as one takes the limit of the large sys- 
tem size L x . Therefore the linear response theory should 
be valid in a sufficiently lar g e s y stem. I now invoke the 
argument used in Refs. |4i ll2lll8j . Namely, as the Hamil- 
tonian is translationally invariant for any value of $, the 
total momentum P x of the system does not change during 
the process. It can be most simply seen in the Heisenberg 
equation of motion dP x /dt = i[H,P x ] = 0. However, at 
t = T, the Hamiltonian H.(<&o) with a unit flux quan- 
tum inserted is not identical to the original one 7i(0), 
although they must have the same spectrum. The large 
gauge transformation 



U 



(8) 



where np is the particle density operator at site r, brings 
the Hamiltonian W^o) back to the original form 7i(0). 
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In the course, the total momentum P x of the system is 
also changed by AP X — 2wvL y x . . . Ld because of 

U^TJJ =T x e iAP % (9) 

where T x = e lPx is the lattice translation operator to x 
direction. The momentum change AP X is different from 
for an incommensurate filling [y not an integer). I note 
that this argument does not require the process to be 
adiabatic. Even for a non-adiabatic time evolution and 
a complicated final state \^>(T)), the momentum is still 
exactly conserved before making the large gauge trans- 
formation, thanks to the translation invariance of the 
Hamiltonian. 

Let me denote E^p x be the energy of the lowest excited 
state with total momentum AP X . Then D is bounded 
from below as D > e 2 L 2 EAp, c /(4:irV). Considering the 
case $(T) — 2irn/e 2 , we can obtain similar bounds on 
EnAPx for n = 1, 2, . . . q - 1. 

This expression suggests the nontrivial dependence of 
the transport property and/or the energy spectrum on 
the aspect ratio |l5|. In taking the thermodynamic limit, 
here I want to keep an isotropic limit. Namely all the 
length Lj are taken to be of same order of magnitude L 
and then L is sent to the infinity. In that limit, we obtain 

D = D> const. L 2 ~ d E n Ap x , (10) 

where d is the dimensionality of the system. This is the 
main result of the present letter. Let me discuss its con- 
sequences for the systems in different dimensions d sepa- 
rately, 
d = 1 

If the system has an incommensurate filling and still 
an insulator (D = 0), 

£nAP x =o(l), (11) 

where n — 1,2, ...q— 1 and o(l/L) means vanishing 
faster than 1/ L as L — > oo. It means that, if there is a gap 
between the ground state(s) and the excited states, there 
are g-fold (nearly) degenerate ground states ^o^) with 
different momenta Pq + 2nkv. (k = 0, 2, . . . , q— 1). Intu- 
itively, the system could acquire a gap by spontaneously 
breaking the translation symmetry to become an insula- 
tor. Then generally we expect g-fold degenerate ground 
states below the gap. This is consistent with the above 
bound. On the other hand, the result l|ll|) does not rule 
out the possibility of a gapless incommensurate insulator. 
Still, the bound (|11J) requires that its finite-size spectrum 
of such a system should not scale (entirely) as 1/L, al- 
though it is the case when the system is described by a 
conformal field theory. 

On the other hand, Resta initiated the argument 0, 
Hif that the phase of the ground-state expectation value 
of eq. © 

z = CW*o>, (12) 



can be interpreted as the polarization which is well de- 
fined only in an insulator. Namely, the system is a con- 
ductor if z = and an insulator if z ^ in the ther- 
modynamic limit. That this approach is rather closely 
related to the LSM argument is pointed out only recently 
by Nakamura and Voit 0] . Because of the LSM argu- 
ment, U\^q) belongs to a different momentum eigenvalue 
than the ground state and thus z = 0. This naively im- 
plies that the system is always a conductor, which is of 
course not true. Observing this, Aligia and Ortiz 
proposed to use z^ = (^o\U q \^o} instead of z, for a 
fractional filling v = p/q. This indeed fits with the 
picture of spontaneously broken translation symmetry, 
which implies the expectation value of the position op- 
erator to be multi-valued. In such a case, the physical 
(pure) ground state in the thermodynamic limit should 

be taken as - Y,k e ^ jhv 1*0*°)- The original or- 

der parameter z can be non-vanishing with respect to 
More analyses would be needed in this direction 
to examine the implications. 

In any case, the present result implies that it would 
be more 'difficult' to make an insulator out of a system 
with highly incommensurate filling (i.e. large q) as there 
are large number of restrictions. Even when the system 
is conducting, requiring the Drude weight to be finite 
we obtain E n &p x = 0(1/ L) if the filling is incommen- 
surate. While this is nontrivial, the standard LSM ar- 
gument 0, 0, 0| already gave this bound irrespective of 
the conduction property. On the other hand, eq. Ijllfl is 
a stronger result which applies only to an insulator. The 
bound (jl 1|> certainly should not apply to a conductor, as 
there are many counterexamples, 
d = 2 

If the system has an incommensurate filling and still 
an insulator (D = 0), 

E nA p x = o(l). (13) 

Namely, they should vanish (at some rate, possibly very 
slowly) in the thermodynamic limit. In particular, when 
there is a finite gap between the ground states and the 
excited states, the ground states must be g-fold degen- 
erate. On the other hand, a gapless insulator phase is 
also consistent with the above bound. This conclusion is 
essentially the same as obtained in Ref. Q assuming that 
the gap does not close due to the AB flux. Here, that 
(questionable) hypothesis is replaced by the assumption 
that the system is an insulator. For a conductor, requir- 
ing the Drude weight to be finite only yields a trivial 
bound E n Ap x = 0(1), Unfortunately, no useful result 
can be obtained for d > 3 with the present argument. 

Straightforward generalization to the multi-species 
cases is possible in parallel to Ref. 0]. For example, let 
me define z/f and v± respectively as up-spin and down- 
spin electron density, and v = + v^. One can consider 
the "real" electric field which couples equally to electrons 
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with either spin, to derive similar constraints if v is frac- 
tional. For the half-filled Hubbard model, this does not 
give any constraint as v = 1 is an integer. However, in 
this case one can still consider a fictitious electric field 
which couples to up-spin electrons only. Then apply- 
ing the same argument as before, one can deduce the 
bounds (II ip> and (| 1 31) respectively for one and two di- 
mensions, if the system is "insulating" with respect to 
the up-spin electrons. For the standard Hubbard model 
in one dimension at half-filling, the exact Bethe Ansatz 
solution shows that all the low- lying states scale as 1/L 
for any value of the Hubbard interaction U. This is actu- 
ally consistent with the above result, because for U > 
the system is insulating with respect to the (real) charge 
but is a "conductor" in terms of spin. For U < 0, it is 
a "spin insulator" but is a charge conductor. Therefore, 
the system does conduct "number of up spin electrons" 
for either signs of U. 

To summarize, in this Letter I have discussed a rela- 
tion among the commensurability, the conduction prop- 
erty of the system, and the finite-size energy spectrum, 
combining the ideas from Refs. 0, an d those from 
Refs. P, 0, 0, 0- As a consequence, I have derived an 
upper bound for energies of some low-lying states for an 
insulator with an incommensurate filling in one and two 
dimensions. It is amusingto observe that several of very 
fundamental concepts 0, LI LLfl LLij LLJ> III! m condensed 
matter physics, some of them put forward about 40 years 
ago, are rather deeply connected with each other, and to 
the simple notion of the commensurability. 

As a final remark, I re-examine the validity of the "gap 
protection" hypothesis made in Ref. |j] . Suppose the sys- 
tem has a finite gap (at zero AB flux) above the ground 
state(s). If the hypothesis is valid (i.e. the gap remains 
finite for any value of the AB flux) , the present argument 
concludes that the system is an insulator (D — 0) . There- 
fore, if there is a conductor (D > 0) with a finite gap, it 
contradicts with the hypothesis. According to Ref. |l5| . 
in the presence of the gap, the superfluid weight D s is 
identical to the Drude weight D, thus the conductor with 
a finite gap automatically means a superconductor. 

In fact, in the presence of the gauge (Coulomb) in- 
teraction, a superconductor generally has a finite mass 
gap, due to the Anderson-Higgs mechanism |20j. Thus, 
the gauge (Coulomb) interaction, being long-ranged, is 
enough to invalidate the hypothesis. If there is no 
long-range Coulomb interaction, a "superconductor" (or 
rather a superfluid) is generally associated to off-diagonal 
long-range or quasi-long-range order, and thus to gapless 
excitations. Hence the hypothesis may still be valid for 
a systems with short-range interactions (and hoppings). 

Recently, Misguich et al. Q verified the "gap pro- 
tection" hypothesis numerically for the RVB spin liquid 
state in d — 2. They actually proposed an even stronger 



hypothesis (at least for the spin liquid state in d = 2) 
that the ground-state energy E(&) is independent of the 
AB flux $ in the thermodynamic limit. This "flatness" 
hypothesis in Ref. [|J may be more naturally understood 
with respect to the conduction property. Namely, if the 
groundstate is always below the gap, its energy E(&) 
must satisfy E(& ) = E(0). Thus, if E(§) varies with <E>, 
it must be convex at one point while it is concave at an- 
other. According to eq. JIJ) the system must then change 
its character from a superconductor to an insulator due 
to the insertion of the AB flux. Excluding such an unnat- 
ural behavior leads to the "flatness" hypothesis in d < 2. 
On the other hand, this argument does not support the 
"flatness" hypothesis for d > 3. 
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